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ABSTRACT 

We consider in more detail a model previously proposed for the hard X-ray (> 
10 keV) emission observed from the supernova remnant Cas A, whereby electrons 
are accelerated by lower-hybrid waves and radiate bremsstrahlung. We consider 
both cold and thermal plasma limits of the modified two-stream instability that 
generates the lower-hybrid waves, and by studying time dependent ionization 
balance for various components of the Cas A ejecta and shocked circumstellar 
medium, find locations within the shell where one or other of the instabilities 
may occur. Either instability can be effective, with the cold plasma instability 
imposing fewer constraints on the shocked reflected ion population responsible for 
exciting the waves. The instability must be located in the ejecta shocked at the 
earliest times and therefore closest to the contact discontinuity where magnetic 
fields are expected to be the strongest. The energy deposited in this ejecta by 
collisions between accelerated and ambient electrons is broadly consistent with 
that required to reheat this ejecta to the observed temperature of ~ 4 x 10^ K. 

Subject headings: acceleration of particles — radiation mechanisms: non-thermal — 
shock waves — supernova remnants 



1. Introduction 

The supernova remnant Cassiopeia A has long intrigued astonomers across the entire 
electromagnetic spectrum. Long known as the brightest radio source in the northern sky 
(Reber 1944; Ryle & Smith 1948), it has also recently been detected in TeV 7 rays (Aharonian 
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et al. 2001), and is most likely the remnant of a putative supernova observed as a 6tli 
magnitude star by Flamsteed in 1680. The discovery of hard X-ray continuum emission 
extending out to 100 keV has stimulated interest in electron acceleration in its shock waves. 
As discussed in a recent paper (Laming 2001), the original interpretation of this emission as 
synchrotron radiation from highly relativistic electrons (Allen et al. 1997) deserves further 
study Unlike in the cases of SN 1006 (Koyama et al. 1995), G347.3-0.5 (Slane et al. 1999; 
Koyama et al. 1997), and now also G266.1-1.2 (Slane et al. 2000), where the X-ray continuum 
is strong enough to completely mask the line emission that should also be present, the 
continuum in Gas A is merely a high energy "tail" to the otherwise Maxwellian thermal 
bremsstrahlung continuum, as shown in Figure 1. (Figure 1 also shows some models discussed 
in detail below). It is thus possible that such emission could be the result of electrons 
accelerated by coUisionless processes other than first order Fermi acceleration at shock fronts. 

Such a scenario is possible because the shock transition occurs on a length scale much 
shorter than particle mean free paths to Goulomb collisions. Assuming that the gas pressure 
>> magnetic pressure in the downstream plasma, the temperature Ta reached by particle 
species a with mass is given by conservation of energy, momentum and particles across 
the shock discontinuity by 

3 

^B^a = Y^^o.^s (1) 

for shock velocity Vg- Thus in the limit of true coUisionless plasma, the particles will have 
temperatures proportional to their masses, i.e. protons will have a temperature mj/me = 
1836 times that of the electrons. Obviously in realistic conditions Goulomb equilibration will 
take place, with rate (Spitzer 1978, see also Appendix B) 

. -0.13ZV.^ (2) 

dt at!'^ 

in c.g.s. units, where AT — — T^, and the ions have charge Z and atomic mass A. 
However coUisionless processes involving the generation of plasma waves at the shock front 
and the associated electron heating by Landau damping or stochastic processes can operate 
on a much faster timescale. Recent simulations of such processes are given by Dieckmann 
et al. (2000) and Shimada & Hoshino (2000) for quasi perpendicular shocks (magnetic field 
in the plane of the shock front) and by Bykov & Uvarov (1999) for a quasi parallel shock 
(magnetic field parallel to shock velocity vector), where shock heated electron distributions 
characterized by Maxwellians at low energies going over into power law distributions at high 
energies are frequently seen. 

Although these results appear qualitatively consistent with the observed X-ray spectra 
of Gas A, the situation is complicated by the fact that optical (Ghavamian et al. 2001) 
and UV (Laming et al. 1996) diagnostics for coUisionless electron heating at fast shocks 
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( admittedly in other supernova remnants, Tycho and SN 1006) give essentially no electron 
heating beyond that implied by the shock jump conditions, equation (1). Similar results have 
been found in X-ray spectra for E0102.2-7219 (Hughes et al. 2000), SN 1987A (Burrows et 
al. 2000), and SN 1993 J (Fransson, Lundqvist, & Chevalier 1996). Laming (1998, 2000, 
2001) speculated that the extrapolation involved from the modest Mach numbers (of order 
1-10) of the simulated shocks to the very high (> 100) Mach numbers for forward and 
reverse shocks in supernova remnants might not be valid. The reasoning behind this is that 
the higher Mach number shocks become progressively more turbulent (c.f. Tidman & Krall 
1971; Tokar et al. 1986; Cargill & Papadopoulos 1988). Although this appears to have no 
effect on the simulations, these have only one spatial dimension. In two or three dimensions 
the conjecture is that the shock reflected ions, while still quasi-monoenergetic, no longer 
behave as a "beam" , now having a much larger distribution in pitch angle. The consequence 
of this is that the preshock instabilities heating the electrons are now no longer reactive in 
nature but kinetic, and have growth rates low enough that insignificant wave growth and 
hence electron heating can occur before the preshock plasma is overrun by the shock front 
itself. Reactive and kinetic instabilities are discussed more fully in Appendix A.l. 

For this reason Laming (2001) investigated a slightly different mechanism; the excitation 
of lower hybrid waves and their consequent acceleration of electrons as a result of secondary 
shocks passing through the presumed region of strong magnetic field in Cas A. Lower hybrid 
waves are electrostatic ion oscillations, which can occur in magnetic fields strong enough that 
the electron gyroradius is smaller than the lower hybrid wavelength. Due to this necessary 
magnetization of the electrons, the waves propagate preferentially across magnetic field lines. 
The parallel component of the wavevector k\\/k < ujpi/ujpe. Since uj/k± ^ (me/rrii)^^'^ uj/k\\ 
the wave can simultaneously be in resonance with ions moving across the magnetic field and 
electrons moving along magnetic field lines, which facilitates collisionless energy exchange 
between ions and electrons on timescales much faster than that associated with Coulomb 
collisions. A cold plasma theory for lower-hybrid waves is given in the Appendix of Laming 
(2001). They have also been previously discussed in connection with cometary X-ray emis- 
sion, supernova remnant shock waves and Advection Dominated Accretion Flows (ADAFS) 
(Vaisberg et al. 1983; Krasnosel'skikh et al. 1985; Bingham et al. 1997; Shapiro et al. 1999; 
Bingham et al. 2000; Begelman & Chiueh 1988), and observed in situ together with acceler- 
ated electrons at Halley's comet (Gringauz et al. 1986; Klimov et al. 1986). 

Lower hybrid waves are especially attractive to consider due to the strong magnetic 
fields inferred to exist in Cas A (e.g. Longair 1994), near the contact discontinuity between 
ambient circumstellar material shocked by the blast wave and supernova ejecta heated by the 
reverse shock. The situation is shown schematically in Figure 2. The explosion drives a blast 
wave into the surrounding circumstellar medium. As this begins to sweep up more plasma 
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behind it and decelerates, a reverse shock develops which moves back into the ejecta (in a 
Lagrangian sense; it is still expanding in Eulerian coordinates), with a contact discontinuity 
between the reverse shocked ejecta and forward shocked circumstellar medium. The dense 
shocked ejecta being decelerated by the much less dense shocked circumstellar medium is 
subject to a Rayleigh- Taylor instability, which is the likely mechanism for the magnetic field 
amplification. Also shown in Figure 2 are two of the various classes of knots or clumps of 
material in Cas A that are particularly conspicuous in optical observations. Quasi-stationary 
flocuUi (QSFs) are seen ahead of the forward shock and fast moving knots (FMKs) ahead of 
the reverse shock. As these shocks run into the density contrasts represented by these struc- 
tures, they split into a transmitted shock that propagates through the clump and eventually 
destroys it, and a reflected shock that moves back towards the contact discontinuity. It is 
these reflected shocks that in our model generate the lower hybrid waves and accelerated 
electrons in the high magnetic field region. 

In this paper wc extend the analysis of Laming (2001) to more realistic thermal plasma 
conditions in Cas A. These are taken from time dependent ionization balance calculations 
based on analytic expressions for the hydrodynamics from McKee & Truelove (1995), and 
updated in Truelove & McKee (1999). Section 2 develops the plasma theory for lower hybrid 
waves in thermal conditions. Section 3 describes the ionization balance calculations with and 
without the extra coUisionless electron heating. Since the mechanism depends on having 
strong magnetic fields in the shell of Cas A, amphfied presumably from seed fields by a 
Rayleigh- Taylor instability, the history of which over the lifetime of the remnant is presently 
unknown, it is not possible to calculate definitively the hard X-ray spectrum. However if 
this model is basically correct, a number of inferences can be made about the nature of the 
plasma in Cas A, and these are discussed in section 4, followed by conclusions in section 5. 
Many of the mathematical details in section 2 and 3 are given in the Appendices. 



2. Cas A and Lower Hybrid Waves Revisited 

In a previous paper (Laming 2001) a theory of the growth of lower-hybrid waves was 
described in the cold plasma approximation (i.e. where the wave frequency a; >> \/2kvi 
and the electron cyclotron frequency fig >> k_i_Vf.\ v^. and Vi are the electron and ion thermal 
velocities respectively), where ions reflected from shocks propagating throughout the shell of 
Cas A generate waves through a modiflcd two stream instability at a perpendicular shock. 
In this case, waves moving at angle 6 to the magnetic field direction given by cos 6' = uipi/ujpe 
are preferentially excited, where ujpi and uj^^ are the ion and electron plasma frequencies 
respectively. This is because waves at this angle have a group velocity away from the shock 
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equal to the shock velocity itself (McClements et al. 1997), and hence stay in contact with 
the shock reflected ions for arbitrary lengths of time. This allows large wave amplitudes 
to develop, even though the growth rates may be small. These waves are excited by ions 
returning to the shock front with velocity 2vs directed along to shock velocity vector. As 
the shock becomes quasi-parallel this feature disappears, since the wave group velocity is 
no longer in the same direction as the shock velocity. The generalization of this theory for 
thermal electrons and ions is given in Appendix A. Where a; >> \/2kvi, it is found that 
outgoing ions can now excite the necessary waves, propagating at direction cosines essentially 
the same as those found in the cold plasma limit for the returning ions, i.e. x'^ ~ cUpJcUp^. 
In the opposite case, a; << y/2kvi, again only returning ions can excite the waves, but for a 
given wavevector do so at a variety of small direction cosines. 

In the hot ion limit (a; << y/2kvi) waves can be generated with a range of direction 
cosines and frequencies, so we now proceed to determine the wave direction cosine and 
frequency at which the net growth rate is maximized in this limit. The growth rate due to a 
Maxwellian ion distribution moving with bulk velocity U in the limit a; << \/2kVi is given 
by a generalization of equation (A12) of Laming (2001), which we reproduce here. 
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where / = ^ /q^"" 4 (^) cxp (^^) v^dv^, vl = k^Tjm,, = k^T,/m,, Jq (^) 



is the Bessel function of the first kind of order zero and Qe is the electron gyrofrequency. 
The wave frequency in this limit is given by equation (A7) with (p = uo'^/k'^vf + ... 

(4) 



l+u;lJk^vf + {l-I)u;^Jk^vr 
The growth rate must exceed the Landau damping rate, Jld, of the wave due to thermal 

— * — * 

ions, given by equation (3) with k ■ U — 0. 

Since a; << V2kvi, we may put k-U » uj and find uj for which ^ — ^ld is maximized. 
Substituting for k\\/k in terms of cu in the expression for 7 — ^yLD, we find 
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Setting this equal to zero we find maxima or minima of uj given by the cubic equation where 
k-U = —kvs/a, with Vg the shock velocity and a = 1/ (2 cos /3) , where /3 is the angle between 
dujjdk^ and U ; 
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+ OJ 



- 2k''vf^ J exp {-W/a^) - 2k''vf)]^ - ^k^vf^^ exp {-M^a^) = (6) 



where M = Vg/ (\/2?;i), the Mach number divided by \/2. For —kvg/a >> kvi this has two 
solutions at u; = \/2kvi corresponding to minima and a maximum at u; = 0. Looking for 
more exact solutions near a; = we drop the term in uj^ and solve the resulting quadratic 
equation to get 



u = -a—V2kvi- '-— — 
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The value of x is then given by 
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and the growth rate by 
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It is shown in the Appendix that under a broad range of conditions, —l<a< —1/2. 
Substituting equation (7) into equation(3) is it easy to see that the growth rate is maximized 
as q; — > —1, and that waves grow fastest for the maximum allowable k. In realistic conditions 
kmax ~ ^e/'^e, siuce the electron gyroradius must be shorter than the wavelength for lower 
hybrid waves. Then for M = 1 the growth rate is ~ 0.04 {n'^cUp^ / riiUJp^^ (QeVe/vi) I / (1 — /) 
rad . For M = 2 the growth rate is 2 orders of magnitude lower. Thus for n[ ~ 0.2nj 
the growth rates are of order 10^^ — 10"*^ x {ve/vi) I / (1 — /) rad s^^ for plasma parameters 
similar to those in Figure 8 [tOpe = 1.8 x 10^, ujpi = ?>x 10^, VL^ = 1.8 x 10^ corresponding to 
B = 1 mG). 



The accelerated electron distribution function is given by the generalization of equation 
(A16) in Laming (2001) 



{—n[kvs exp (— M^/ct^) /a — riiUj) OJ. 
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where constant of integration, the maximum electron velocity where /g {Vm) = 0. 

Integrating over ?,'[[ in the range Vm, and making the approximations cu « ^/2kvi, 
k ■ U = kvs the accelerated electron density is 



^ 8M 



2 4 2 



- 7- 



For reference, the accelerated electron density in the opposite limit, a; >> \/2kvi is 
given by (Laming 2001, equation (A17)) 



with X ~ LOpi/uOpe. Taking uo to be the lower hybrid frequency, equation (11) is valid for 
Ti » 2.3 {B/lmGf x lOVue K and equation (12) for « 2.3 {B/lmGf x lOV^e K. In 
between these limits Landau damping by thermal ions is sufficiently strong to quench the 
instability. 

To summarize, we have calculated wave dispersion relations and kinetic growth rates in 
each of the limits u >> \fikvi and uj « \/2kvi. Since these growth rates are in general 
lower than the ion cyclotron frequency, fli, we also calculate the group velocity of the waves, 
and by setting this equal to the shock velocity find the angle to the magnetic field at which 
propagating waves can remain in contact with the shock refiected ions and hence grow to 
large amplitudes over essentially arbitrary periods of time. For these waves we calculate the 
growth rates and hence densities of accelerated electrons. 

We have revised the velocity range over accelerated electrons may exist from —v^ — ^ Vm 
in Laming (2001) to — > due to the following consideration. Karney (1978) has derived 
a value for the maximum electric field above which the wave becomes stochastic and ion 
heating also commences, given hy E — B {VLi/uj)^^^ uj / 'i/k±_ — B (Qi/uj)^^^ Vs/2. The kinetic 
energy acquired by an electron in time t ~ (the time available for electron acceleration 

in the precursor before the shock hits) is (mf/2me) {E / B)^ / Z"^, which for relevant values of 
X and uo evaluates to energies less than 1 keV for Cas A. This is valid for a true perpendicular 
shock. However it is easy to see that an electron traveling along a magnetic field line inclined 
at a small angle to the shock front (a quasi-perpendicular shock) can easily stay ahead of 
the shock and be accelerated to large energies, but may do so traveling in one direction only. 

The calculation of the accelerated electron distribution function assumes that the elec- 
trons are accelerated purely by Landau damping of the lower hybrid wave, which is likely 
to be the case for the non-relativistic electrons envisaged here. However since lower hybrid 
waves exist with u}/k\\ ~ c, electrons can in principle be accelerated to relativistic energies, 
if they can be trapped in the wave generation region for long enough. McClements et al. 
(1997) argue that this is indeed the case for quite a wide range of shock Alfven Mach num- 
bers, based on the intensity at which lower hybrid waves saturate in simulations and the 
resulting value of the momentum diffusion coefficient for electrons, and that lower-hybrid 
waves may provide an electron injection mechanism for subsequent Fermi shock acceleration 
to cosmic ray energies. 




(12) 
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3. Reverse Shock Ionization Structure 

3.1. Initial Conditions and Basic Equations 

We model the hydrodynamics for the Cas A shell assuming uniform density ejecta and 
ambient medium. We use expressions given by McKee & Truelove (1995) for the reverse shock 
radius and velocity with respect to the expanding ejecta, Rr and Vr, which connect between 
the early ejecta dominated, and late Sedov- Taylor limiting behavior. The generalization 
of such models to non-uniform densities is given in Truelove & McKee (1999). We choose 
parameters for the hydrodynamic models of Cas A as follows. The blast wave velocity is 
inferred to be 5200 km s~^ (Vink et al. 1998; Koralesky et al. 1998). The expansion factor 
{dlogRh/dlogt where Rb is the radius of the blast wave) found by these authors to be 
Tj — 0.64 suggests that the remnant is slightly more than midway through the transition 
between the ejecta dominated and Sedov- Taylor dominated phases, where values of 77 of 1.0 
and 0.4 respectively would be expected. Hence we expect Rst — 2.23 (Me/Mo)^*^^ rig ^''^ 
pc and tsT = 209£'5/''^ {Mg/MQ^^^ n^^^^ years, the transition blast wave radius and time 
(McKee & Truelove 1995), to be similar to or perhaps slightly smaller than the current values, 
2.64 pc and 320 years respectively. In these relations no is the preshock hydrogen density, 
-E51 is the explosion energy in units of 10^^ ergs, and M(,/Mq is the ejecta mass in units 
of the solar mass. The ionization age for this shock front is determined to be in the range 
10^^ - 2 X 10" cm-=^s, (Vink, Kaastra, & Bleeker 1996; Favata et al. 1997), which combined 
with the age of 320 years gives a preshock electron density of 2.5 - 5.0 cm~^, assuming a 
shock compression of a factor of 4. Taking Rst — Rb — 2.64 pc gives M^/Mq in the range 
4-8. Assuming tsr — 320 years gives £'51 in the range 2.4-4.8. Prom ASCA observations, 
Vink, Kaastra, & Bleeker (1996) find Me/Mo ~ 4. Favata et al. (1997) analyzing BeppoSAX 
spectra find 2 < M^/Mq < 4, with the uncertainty stemming from different assumptions 
about the role of the non-thermal continuum. The small ejecta mass is consistent with 
extensive mass loss as stellar winds by the progenitor. Both references find O to be the 
dominant element composing the ejecta, with rather little of the expected products of O 
burning evident, suggesting further that these are confined in a compact stellar remnant. 
The emission measure quoted by Laming (2001) corresponds to AMq of O (if Ue — 10 cm~^), 
so in the following we will take £^51 = 2, M^/Mq — 4, and Uq — 3 cm~^. This values are 
similar to those considered by Mochizuki et al. (1999), with the exception of a smaller hq. 
It appears that they chose uq consistent with the modeling of Borkowski et al. (1996), who 
modeled the deceleration of the blast wave upon encountering a shell of circumstellar material 
postulated to be at the interface of red and blue supergiant presupernova winds, in order 
to match the measured radio expansion of ~ 2000 km s~^. The more rapid expansion seen 
in X-rays (Vink et al. 1998; Koralesky et al. 1998), generally assumed to be more directly 
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associated with the actual motion of the plasma, suggests a lower circumstellar density as 
used here. Borkowski et al. (1996) do demonstrate in their simulations one important feature 
relevant in this work, namely the fact that the blast wave (and presumably also the reverse 
shock) split into transmitted and reflected shocks upon encountering density contrasts, and 
that these secondary shocks will propagate back across the presumed region of high magnetic 
field near the contact discontinuity. 

The blast wave of Cas A is observed at a radius of 160", while the bright ring of X-ray 
emission is at a radius of 110" (Vink et al. 1998; Gotthelf et al. 2001). This suggests that the 
bright ring is in fact emission from the ejecta heated by the reverse shock, since according 
to the self-similar models of Chevalier (1982), for a blast wave radius of 160", the radius of 
the contact discontinuity should be in the range 120" -140". Behind the reverse shock, the 
density Ug of ions with charge q is given by 

— {Cicyn,q-inq^i — Cion,q'>T'q) + [{Crr,q+1 + C'dr,q+l) '>T'q+l ~ {Crr,g + C'd^.q) ^q] (13) 

where Cion,q, Crr,q, Cdr,q are the rates for electron impact ionization, radiative recombination 
and dielectronic recombination respectively, out of the charge state q. These rates are the 
same as those used in the recent ionization balance calculations of Mazzotta et al. (1998), 
using subroutines kindly supphed by Dr P. Mazzotta (private communication 2000). The 
electron density Ue is determined from the condition that the plasma be electrically neutral. 
The ion and electron temperatures, Tj and are given by generalizations of equation (2) 
for the case when many charge states of a particular element are present (see Appendix B) 

and 

dT,^ (r,-re)E,gX/(g+l) T, fdriA _^_dQ 

dt at!'' T^qUq V dt ZnM df ^ ^> 

Here A is the atomic mass of the ions in the plasma. The last term dQ/dT represents plasma 
energy losses due to ionization and radiation. Radiation losses are taken from Summers & 
McWhirter (1978). These authors tabulate the power emitted by each ion of the elements 
C, O, Si, Ar, Fe and Mo as a function of temperature, allowing an exact treatment of 
the ionization non-equilibrium aspect of the power losses for these elements. Losses for 
other elements can be readily obtained by interpolation along isoelectronic sequences. At 
each time step Tj and are modified by a further factor exp (— 4vea;/3r) and rig and the 
Uq by exp {—2vex/r), coming from the adiabatic expansion of a spherical shell of plasma 
with volume V — Airr^dr where dr is held constant as r increases. In this approximation. 
Vex, the expansion velocity, which is constant for all ejecta, is given by the lesser of Vgx — 
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R^/tr — 0.75iv and v^x — 0.75i'6, with tr being the reverse shock onset time, and being the 
blast wave velocity. This is a reasonable approximation for the transition period between 
ejecta dominated and Sedov- Taylor evolution in a supernova remnant (sec the "hollow" 
blast wave approximation in Ostrikcr & McKcc 1988). More detailed work could use the 
self similar solutions appropriate to the various regions and epochs of the remnant evolution 
(e.g. Chevalier 1982; Hamilton & Sarazin 1984), but even then the inevitable Rayleigh- 
Taylor instabilities are still difficult to treat. The radius r of the plasma element is given 
hy r — Rr -\- Vexdt, and we neglect the fact that the rear side of Cas A is around 4pc 
further away from us than the front side, and thus is observed about 13 years further back 
in time. Equations (14) and (15) are integrated forwards in time from initial conditions 
given by the shock jump conditions, equation (1), for the electron and ion temperatures 
and by the hydrodynamic model discussed above for the densities. All ions are assumed to 
be singly charged at the time of shock passage, and all ions of the same element arc taken 
to have the same temperature. Our results are insensitive to the exact initial charge state 
distribution since we arc mainly interested in what happens many years after the initial 
postshock ionization has ceased. In comparisons with the calculations of Mochizuki et al. 
(1999) for pure Fe ejecta we found agreement to ±10% or better for time dependence of the 
electron temperature over the lifetime of Cas A. 



3.2. Ionization Structure with no Collisionless Electron Heating 

Figure 3 shows the evolution of the ion temperature (dashed lines) and the electron 
temperature (solid lines) for O ejecta heated by reverse shock encounters at various times 
in the evolution of Cas A. The highest electron temperature, Tg — 2.7 x 10^ K is found in 
plasma shocked about 150 years after the supernova. Table 1 gives a summary of the present 
day values of plasma parameters for different shells of the ejecta shocked at different times 
after the initial explosion. In Table 2, and illustrated in Figure 4 are similar temperature 
plots but for circumstellar plasma with elemental composition dominated by N, He and 
H (in mass ratios 0.02:0.49:0.49 respectively) shocked by the blast wave. Here the highest 
electron temperature is 7.4 x 10^ K, again for plasma shocked 150 years after explosion. This 
value for the electron temperature is consistent with that actually observed directly behind 
the blast wave (U. Hwang 2001, private communication), but higher than that observed 
in spatially integrated spectra. However we believe that the reverse shocked ejecta is the 
dominant contributor to the observed X-ray emission. Besides the morphological arguments 
given in the previous section, the blast wave as modeled does not provide enough emission 
measure to account for the observed emission by a factor of at least 4, remembering that 
this plasma being composed of less highly charged ions radiates less bremsstrahlung per unit 
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emission measure. Although one could increase the circumstellar density and ejecta mass in 
the model to match the observations, the explosion energy would also have to be increased 
as well to to maintain the consistency with the observed blast wave velocity and radius. A 
problem might still remain however with the likely overproduction of O Vlll emission from 
the dense ejecta cooling at temperatures around 10® K, compared with observations (e.g. 
Bleeker et al. 2001). Only ejecta which encountered the reverse shock early in the evolution 
of Gas A (i.e. at times less than 100 years after explosion) has sufficient emission measure 
to account for the observations. 

High values of the ion temperatures in our simulation result simply from the choice 
of the shock jump conditions, equation (1), as our initial conditions. Some observational 
support for very high postshock ion temperature exists. Ghavamian et al. (2001) sees broad 
H a. profiles in optical spectra of Tycho, RGW 86 and the Gygnus Loop, commensurate with 
the shock velocities derived in these remnants. Further, Raymond, Blair, & Long (1995) 
observing the NW limb of SN 1006 in the UV see broad line profiles for G IV, N V, and O 
VI, indicating that the heavy ions have similar velocity distributions to each other and to 
the protons, and temperatures proportional to their masses consistent with equation (1). 

So we proceed to place the dominant X-ray emission and the electron acceleration in 
the ejecta. One immediate, but easily solved problem is that according to the ionization 
balance calculations the temperature of the ejecta shocked at early times should be below 
10® K, not ~ 4 X lO'^ K as observed. However, if the hard X-rays arc due to bremsstrahhmg 
from suprathermal electrons, one should expect just such a temperature deficit in order to 
be able to accommodate the inevitable Coulomb heating of the ambient plasma due to the 
accelerated electrons. Gas A emits ~ 10^^ ergs s~^ above 10 keV. The ratio of radiated 
bremsstrahlung power to that lost to Goulomb collisions is given approximately by (see e.g. 
Katz 1987) ZgP^afJ (Sv^logA) where ~ 1 is the bremsstrahlung Gaunt factor, (3 — v/c, 
and afs is the fine structure constant. For Gas A, Z = 8, /3 ~ 0.2, and log A ~ 40 giving 
a ratio of ~ 10~^. The nonthermal emission above 10 keV thus requires an energy input 
of ~ 10"^° ergs s~^, or a total energy input integrated over the lifetime of the remnant of 
~ 10^*^ ergs. The electron internal energy in an emission measure of UeUoV = 3 x 10^^ cm~^^ 
for no ~ 2 — 5 cm~^ (corresponding to ejecta shocked 50 and 100 years after explosion) 
is 2.bneVk-B,T ~ 10^^ ergs for a temperature of 4 x 10^ K. However the energy required to 
maintain this quantity of plasma at such temperature against energy losses from ionization, 
radiation and adiabatic expansion since the initial reverse shock encounter is more than this, 
~ 10^° ergs, similar in order of magnitude to that available for heating from the accelerated 
electrons. 



3.3. The Inclusion of Accelerated Electrons 



To investigate further the effect of the Coulomb heating of the ambient electrons by 
the accelerated population, we plot in Figure 5 the time dependence of the electron and 
ion temperatures behind the reverse shock under various heating conditions. We assume a 
constant heat input per unit volume fixed at the present day value in ejecta reverse shocked 
about 50 years after the explosion, ~ 4 x 10"^^ ergs s~^cm~^, (see Appendix B.2 with 
He — 50 and n'^ — 2). This heating is taken to "switch on" at different times after reverse 
shock passage. These are taken to be 150, 200, 250, and 300 years. A physical mechanism 
for such a "switch on" is suggested by the morphology of SN 1987A. The blast wave must 
encounter a strong density gradient before reflected shocks can travel back towards the 
contact discontinuity. Thus a certain amount of time may elapse while the blast waves reaches 
this density gradient and the reflected shock propagates back. In SN 1987A the relevant 
density gradient would be provided by the inner edge of the presumed equatorial disk, the 
recombination radiation from which currently produces the inner ring. In Cas A no real 
evidence for such an equatorial disk exists, but shock reflection from a distribution of quasi- 
stationary fiocculi is the main possibility. ^Prom Figure 5 it is apparent that electron heating 
beginning 200-250 years after reverse shock passage (or after 250-300 years in the evolution 
of the supernova remnant, i.e. between curves d and e) will reproduce the current observed 
electron temperature of ~ 4x 10^ K. Figure 1 illustrates this with spectra from the BeppoSAX 
MECS (Boella et al. 1997) and PDS (Frontera et al. 1997) instruments acquired on 1997 
November 26, plotted against a pure thermal bremsstrahlung spectrum at 4 x 10^ K and a 
similar spectrum including an extra 4% of electrons in an accelerated distribution function. 
Also shown are seven data points from a rocket flight on 1968 December 5 (Gorenstein, 
Kellogg, & Gursky 1970). These are shown in more detail in Figure 6 plotted against pure 
thermal bremsstrahlung spectra corresponding to temperatures of 2 x 10^, 3 x 10^ and 4 x 10^ 
K. From Figure 5 it is clear that on the basis of the hydrodynamic models used in these 
calculations, the ambient electron temperature should have increased between these two 
observations by ~ 2 x 10'' K. Depending on how one treats line emission and interstellar 
absorption (both omitted from the theoretical spectra in Figure 5), such a temperature 
increase is marginally consistent with the observations. The ionization age Uet in such a case 
would be ~ 10^^ cm~^s"^, a factor of ~ 2 lower than observations (Vink, Kaastra, & Bleeker 
1996; Favata et al. 1997). Curves b and c in Figure 5, representing electron acceleration 
starting 150 or 200 years after explosion (100 or 150 years after reverse shock passage) give 
values of ngt more consistent with observations, but temperatures in the range 5-6 keV. 
The heating mechanism envisaged here only occurs at perpendicular shocks, and so not all 
the ejecta will be heated. Consequently a thermal conduction energy loss from the heated 
regions might be expected, which would reduce the peak temperature and reduce the degree 
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of temperature increase that might be expected between the observations of Gorenstein, 
Kellogg, & Gursky (1970) and the present day. One interesting feature of curves b and c 
is that the initial period of heating from 10^ K to ~ 2 x 10^ K is very fast, but then slows 
abrubptly. The "plateau" at this temperature is due to the strong peak in the O radiative 
loss curve at this temperature, and the plasma temperature only starts to rise again once it 
has expanded sufficiently to reduce the density so that the coUisionless heating can overcome 
the strong radiative losses. 

We will now investigate in more detail exactly where in the Gas A ejecta a two stream 
instability that provides lower hybrid wave electron heating may occur. In the picture that 
we develop, the X-ray emitting shell of Gas A is filled with shocks propagating in all directions 
as a result of forward or reverse shock interactions with density inhomogeneities. In the hot 
ion limit, where a; << \/2kvi the wave growth rates fall dramatically with increasing Mach 
number, so we can reasonably assume that only shocks with the minimum Mach number to 
be supercritical are important. For plasma (3 (ratio of gas pressure to magnetic pressure) 
of 1 this Mach number is 1.7 for perpendicular shocks, falling to 1.35 for jS = 2 (Edmiston 
& Kennel 1984). A further simplifying assumption will be to assume that the accelerated 
electrons heat the ambient plasma to the observed temperature, 4 x lO'' K, and then that 
the instability continues in marginal stability. This is a lesser criterion on the electron 
energization process than the ambient electron heating discussed above. In such conditions, 
the values of the parameters / and {k±/2) dl/dk± in equation (A7) are ~ 0.5 and ~ —0.25 
respectively. The accelerated electron density is maximized for a given k at —M/a ~ 0.7 
from equation (11). This is a lower value than is realizeable in practice, the minimum value 
being —M/a ~ 1 for high /3 plasmas and short wavelength waves. At marginal stability, 
the electron plasma (3e is given by I3e = SnriekTe/B^ = 2uoiyjnl/c^ = 2. In heavy element 
plasma with highly charged ions the electrons dominte the total plasma f3, so M ~ 1. For 
k = ijjpejc = Qf./ve, —a = 0.5 — 0.6 (from Figure 7), so —M/a ~ 2. Decreasing /?, allowing 
larger k and hence larger (i.e. more negative) a also requires larger M, so —M/a ~ 2 
holds over the range of reasonable values of /3, k, and a. Assuming n'Jui — 0.2 for the shock 
reflected ions, which is comparable to that observed in situ at solar wind shocks (e.g. Gedalin 
1996), the fraction of electrons accelerated is given by n'Jue — 0.05/ vf where Vm — i-^x 10^^ 
cm s""*^ has been used in equation (11). This corresponds to 70 atomic units found as the 
best match to BeppoSAX data in Laming (2001). 

The density of accelerated electrons may be increased over that given by equation (11) 
if the energetic electrons are pitch angle scattered away from the magnetic held direction. 
This scattering mainly results from collisions with ions, while energy loss is due to collisions 
with ambient electrons. The isotropic density of accelerated electrons, n'^iso is given by 
^e,iso ~ f^'e^e,e/'Te,i ^ 4ng. The clectron-elcctron and electron-ion scattering cross sections 
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have the same dependence on the energy of the accelerated electron, and so the isotropic 
accelerated electrons will have the same energy distribution as those directed along the 
magnetic field. This pitch angle scattering boosts the accelerated electron density by a 
factor ~ 5 for the same refiected ion density. Hence an accelerated electron fraction of 4% 
may be possible in regions of the Cas A ejecta where Vi < 2.5 x 10^, i.e. Tj < 6 x 10^ K. 

^Prom Table 2 we see that ejecta reverse shocked up to ~ 100 years after the initial 
explosion may support sufficient electron acceleration to account for the observations. In 
this context it is important to realize, however, that ejecta shocked within 100 years of 
the explosion can have a low enough ion temperature to justify using equation (12) for 
the accelerated electron density, which would increase the accelerated electron fraction by 
around an order of magnitude over that given by equation (11). In general cool plasma 
instabilities are stronger than those in hot plasma, and so our conclusion that electrons can 
be accelerated in ejecta shocked early in the evolution of Cas A will be strengthened. 

Hence our model for the electron acceleration to form the hard X-ray tail by non- 
thermal bremsstrahlung and to provide heating of the ambient ejecta plasma to ~ 4 x 10^ 
K by Coulomb coUisions requires the existence of shocked ejecta that has cooled to give 
ion temperatures well below ~ 10^ K. This requires that secondary shocks passing through 
this ejecta have velocities below ~ 500 km s~^ so as to avoid reheating the ions up to 
temperatures where the instability weakens. This is actually quite fikely, since a shock 
refiected for example from the blast wave moving through the shocked circumstellar medium 
will slow down considerably upon reaching the shell ejecta shell. 

4. Discussion 

We have previously argued (Laming 1998, 2000, 2001) that two stream instabilities 
formed ahead of fast shock fronts in supernova remnants might more realistically be modeled 
as kinetic rather than reactive instabilities. This is observationally motivated, in that the 
degree of electron heating predicted by models based on reactive instabilities (e.g Cargill 
& Papadopoulos 1988) is not observed (Laming et al. 1996; Ghavamian et al. 2001). A 
theoretical justification comes from the fact that high Mach number shocks are turbulent 
and not steady state structures, and reflected ions consequently will not form a "beam" , but 
will inhabit a much wider portion of phase space, giving rise to kinetic instabilities. With 
their lower growth rates, these will generally allow insufficient wave growth ahead of the 
shock before being overrun by the shock itself. In this paper we have avoided this problem 
by focusing on lower hybrid waves formed ahead of the shock that have group velocity away 
from the shock equal to the shock velocity itself, and hence stay in contact with the shock 
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reflected ions for sufficiently long times to allow signiflcant wave growth to occur. 

In this paper we have found that the most efficient electron acceleration occurs for the 
lowest Mach number shocks, just above the first critical Mach number, which is where shock 
is said to become "supercritical" . This denotes the onset of turbulence in the shock struc- 
ture, and ion reflection which becomes progressively more unsteady and "bursty" in nature 
as the Mach number increases (e.g. Tokar et al. 1986). Hence although we considered the 
lower hybrid instability excited by the secondary shocks in the Cas A shell to be kinetic in 
nature since this appears to be a more robust assumption, the results of our calculations 
suggest that in this case reactive instabilities could be quite plausible also. It is also worth 
remarking that the designation "reactive" or "kinetic" comes from the mathematical treat- 
ment of two different limits of the Vlasov equation (see Appendix A), and that physically, 
a continuum of instabilities exists between the kinetic and reactive limits. The further in- 
vestigation of these would require numerical simulations. All instabilities in this paper have 
a minimum wavector kmin ~ ^pe/c arising from electromagnetic effects and a maximum 
kmax ~ ^e/ve from the requirement that the electron gyroradius be smaller than the wave- 
length. In marginal stability kmin ~ kmax allowing us to estimate the magnetic field stength, 
B = m^V(.ujp(./ e = 5.9ujpf. x 10^^, which evaluates to 1.9 mG for ujpe = 4 x 10^ rad s^"*^ in the 
t = 50 year ejecta. ^From Table 1 the ejecta shocked at 50 years after explosion occupies 
a volume ~ 2 x 10^^ cm~^, in order to give the stated emission measure. The magnetic 
energy contained in this volume is 3 x 10^^ ergs, not inconsistent with that deduced from 
the minimum energy synchrotron luminosity argument (Longair 1994). 

These magnetic fields are assumed to be amplified from seed fields by the action of a 
Rayleigh- Taylor instability (though see Lou 1994, for another possibility). Plausible evidence 
of this has recently been demonstrated in imaging spectroscopy of Cas A by Chandra (Hughes 
et al. 2000; Hwang, Holt, & Petre 2000). In the SE quadrant of Cas A in particular, 
emission in Fe L and K shell lines is observed to be outside the emission from the lighter 
element Si, suggesting that the Si burning zone in the progenitor (which produced the Fe) 
has overtaken the O burning zone which produced the Si. This region in the SE is also where 
recent XMM- Newton observations find the hardest continuum (Bleeker et al. 2001), though 
these authors comment that the hard X-ray continuum emission actually is distributed quite 
uniformly over the remnant, which argues in favor of an emission mechanism independent 
of the magnetic field, i.e. bremsstrahlung, and not e.g. synchrotron emission. We have 
previously argued that the secondary shocks that excite the lower hybrid waves arise as 
the blast wave encounters density inhomogeneities in the circumstellar medium producing 
reffected shocks as in Borkowski et al. (1996). However another possibility exists which is 
perhaps more intimately associated with Fc emission and Rayleigh- Taylor instability, viz. 
the reverse shock encountering Fe-Co-Ni bubbles in the ejecta (Li, McCray, & Sunyaev 
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1993; Borkowski et al. 2000). The simulation in the last reference demonstrates that such 
a scenario would produce large amounts of turbulence and vigorous mixing of Fe with the 
overlying ejecta. We also note in passing that the idea of secondary shocks propagating 
through already shocked material, and also being highly magnetic, has been successfully 
applied to ultraviolet spectra of the bowshock in HH 47A by Hartigan et al. (1999). In 
this case the strong magnetic field arises as radiatively cooling gas is compressed by hotter 
material around it. This is also possible but in our view unhkely for the regions of Cas A 
we discuss. 

A common parameter quoted in models for electron heating is the fraction of the shock 
energy that is deposited in the electrons. It is not possible to make such a statement about 
the secondary shocks responsible for the electron acceleration in our model of Cas A, but 
in terms of the blast wave energy an estimate can be made. If the 10^° ergs s~^ required 
to accelerate the electrons comes from secondary shocks that originate from the blast wave, 
this energy is to be compared with the energy ~ AT^R^VsUik-BTi ~ 10^"^ ergs s"-*^ deposited 
by the blast wave in the shocked ions. Neglecting energy going into cosmic ray ions, around 
10% of the blast wave energy ends up in energized electrons to be compared with the ~ 20% 
estimated by (e.g. Cargill & Papadopoulos 1988). 

The electron distribution function is given by balancing the imaginary terms in the di- 
electric response, which give wave growth due to shock reflected ions and Landau damping 
by the accelerated electrons. We argue that secondary shocks maintain this electron dis- 
tribution against the inevitable losses due to Coloumb collisions, and that therefore these 
shocks must pass through a plasma element on timescales significantly shorter than those 
for Coloumb coUisions, which are of order 2-40 years for accelerated electrons of energy 10 
- 70 keV in the ejecta shocked at early times. The accelerated electron distribution will be 
different from that modeled above if processes other than electron Landau damping occur. 
McClements et al. (1997) argue that sufficient lower hybrid wave intensity is likely to ex- 
ist such that electrons may be accelerated by a stochastic mechanism to mildly relativistic 
energies, but that this is more easily satisfied at high Alfven Mach number shocks. 

5. Conclusions 

In Laming (2001) it was shown that in principle a distribution of electrons accelerated by 
lower hybrid waves comprising about 4% of the plasma electrons produced a very good match 
to existing spectra of Cas A from BeppoSAX. This model used cold plasma theory, assuming 
uj » ^/2kvi. In this paper we have extended the plasma theory to treat the opposite limit 
a; << \/2kvi, and by performing calculations of the time dependent ionization balance within 
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the ejecta of Cas A, have identified regions where such wave generation may take place. 

In ejecta shocked ~ 150 years or more since explosion, only the thermal plasma insta- 
bility (a; << y/2kvi) may occur. However this instability appears to be too weak to provide 
significant electron heating to give the hard X-ray emission. In ejecta shocked up to ~ 100 
years since the explosion, either the cold or thermal plasma instabilities (a; >> y/2kvi) can 
occur, depending on the ion temperature. The cold plasma instability is more likely, and is 
sufficiently strong to explain the putative hard X-ray bremsstrahlung. Such electron heating 
and acceleration mechanisms may plausibly occur elsewhere, the main requirement being a 
sufficently strong magnetic field. We have concentrated our efforts on Cas A since it has the 
most conspicuous hard X-ray emission, is well studied in other respects, and is thought to 
have undergone relatively httle elemental mixing (Johnston & Yahil 1984) which renders the 
plasma theory more tractable. 

This work was supported by basic research funds of the Office of Naval Research and has 
also made use of data obtained from the High Energy Astrophysics Science Archive Research 
Center (HEASARC) provided by NASA's Goddard Space Flight Center. 



A. Appendix: Lower Hybrid Waves 

A.l. Reactive and Kinetic Instabilities 

In cold plasma, the dispersion relation for lower hybrid waves, omitting the electromag- 
netic terms is given by (e.g. Vaisberg et al. 1983; McClements et al. 1997; Laming 2001) 

= 1 + ^sin^^ - 4 - '^cos'e^ (Al) 

where uo is the wave frequency, uOp^, and are the electron plasma frequency and gyrofre- 
quency respectively, Upi is the ion plasma frequency and 9 is the angle between the direction 
of wave propagation and the magnetic field vector. As sin^^ 1 this gives uo ^lh 
where ^Ilh = y/^e^i, the geometric mean of the electron and ion gyrofrequencies, known as 
the lower hybrid frequency. In the presence of suprathermal ions the dispersion relation is 
modified by the addition of the following term to the expression for K^; 

+^ f ^-"^^ ^ . ^dh-. (A2) 
where cUp^ f (vi) and n' are the plasma frequency, distribution function and density of the 

— * 

suprathermal ions with velocity Vi, and k is the wavevector with magnitude k. 
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If df /dvi — at uj — k-Vi then the integral may be evaluated in a straightforward man- 
ner. Taking the most extreme limit in this case we put /' = n'5 {vl — vl) for a monoenergetic 

ion beam moving with velocity vl, and the integral evaluates to ujpj (w — k ■ vCj ■ With 
the addition of this term to the expression for the dispersion relation becomes a quartic 
equation analagous to that for the Buneman instability, the complex roots of which (when 
present) give growing or damping waves. Physically the wave growth is driven by the free 
energy in the ion beam as a whole. Such instabilities are known as "reactive" instabilities 
and usually give very fast growth rates. 

In the opposite limit we make the replacement u ^ u + i'j assuming 7 << u; and 
evaluate the integral by going into the complex plane. This procedure is summarized by the 
so-called Landau prescription 

1/ (^tu -k-vlj ^ P (^1/ (^oj - k ■ vi^^ - iTiS (^cu - k ■ . (A3) 

Upon making this replacement, taking imaginary parts and rearranging, an equation for the 
growth rate 7 in terms of the integral f (^k ■ df'/dvj^ 6^ {uj — k ■ ifj^ d^Vi results (see Laming 

2001, for the full expression). In this limit the instability is called "kinetic", and is physically 
due to stimulated Ccrcnkov emission of plasma waves. For wave growth, df /dvi > 0, or in 
other words a population inversion must exist. This makes the analogy with quantum physics 
clear, and indeed one can discuss the various wave-wave and wave-particle interactions in 
terms of the appropriate Feynman diagrams. 

Reactive and kinetic instabilities derive from the different mathematical treatments of 
the term in the dispersion relation accounting for the suprathermal ions. Of course in reality 
a continuum of instabilities exists between these two limits, but further investigation of these 
requires numerical techniques. 

A. 2. The Dispersion Relation in Thermal Plasma 

In finite temperature plasma, the general expression for the longitudinal part of the 
dielectric tensor is (Melrose 1986) 

= 1 + / j::-^ (^) {^ii + / (p) 

= 1 - *^Er-»/r (^)'''-p K«s/2*bt.) rf#4) 
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Here a; is the wave frequency, and k± are the components of the wavevector k parallel and 
perpendicular to the magnetic field B direction, fl^ is the electron cyclotron frequency, v and 
p are the electron velocity and momentum, with subscripts having the same meaning as for 
the wavevector, and Jg is the Bessel function of the first kind of order s. Assuming << 
cu — sfle in denominator, appropriate for lower-hybrid waves, the expression simplifies to 



The term in s = is 



1 pe i 

J- /? iJ? 



2Q2 ( jn^ 



+ k\ 



( ™-e 



X /o^~ (^) exp {-m,vl/2k^T,) 2nv^dvj 



- 



a;2 A;2 k^T^ Jo 



j2 I b^z^ j gxp (-TOe^;l/2/CBT'e) ^nV^dv^, 



(A5) 



(A6) 



while the terms in s 7^ are evaluated by converting the ~^Yli°°- Using the standard 

property of the Bessel functions, I — Jq (z) — 2^^-^ (z), and assuming u; « gives 



a;, 



Tv-i. _ pe 



'|s|>0 



k"^ \knT, 



exp {-mevl/2kBTe) 2Trv±dv±. (A7) 



Including the term for thermal ions, Upjk'^vf ^1 — ? where the usual plasma dis- 

persion function is (p (z) — —z/y/TT exp (— t^) / (t — z) dt, a term unity for the vacuum, 

^ cjpg/ (1 + ujpjk'^c^) , for the electron term 



and inserting an electromagnetic correction, cj^g 
involving k\\ (Begelman & Chiueh 1988; Melrose 1986) the dispersion relation is 

^ < [ik^/e) I (1 + uijc^k^) 
— 5 7, 



where / 



/o^°° -^0 (^n7^) ^^P ( "ct""^ ) = kBTe/rUe and = kBTi/rrii. 



(A8) 



A. 3. Group Velocities 

Proceeding as in the cold plasma case, we differentiate equation (A8) with respect to 
k± to derive an expression for the group velocity; 



dk± 



-a; fcx _|_ ^'^pe(EM) kj_ _|_ ui dl 
k k ~^ k^c^ k ^ 21 dk± 



+ - 



2^pi k± 



dl 
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where = (l + ^pe/c'^')- 

In the limit a; >> \/2kvi, (p / uj'^+?>k*vf / uj*+.... We set duj/dk±_ = Vs = auj/k±_ 

where Vg is the shock velocity and a—1/ {2 cos (3) where (3 is the angle between duj/dk±_ and 

— * 

U, the bulk velocity of the shock reflected ion distribution. With x — k\\/k to lowest order 
in ulJul^, 



A T / i2 T 



+ 



+ 



X^ 



'^pe(EM) 2 afcx "'^pe(EM)-' fe2c2 l+a;2^(l-/)/ifc2^2 

+^pi - (^^vi + i+^2^(i_/)/fc2^2 U. (AlU) 

This looks like a quadratic equation in x^, but we must remember that k^ = k\/ (1 — x^), 
ujpe[EM) ^pe and rewrite in terms of k±, which gives a quartic equation for x'^ . However 
we are looking for solutions x"^ « 1, for which both quartic and quadratic equations give 
the same result to leading order in ^pil^pe{EMy -^^ limit k±Ve/fle 0, 

where 7 ~ 1 — k'j_vl/Ql and {k±/2) dl/dk± ~ 7 — 1 this has solution 

7—^ TT^ Wi/^i(EM)) ■ (All) 

In the cold plasma approximation, '^pe(_BM)/^^'^^ ^ ^^'^ ~ ^pi/^pe ^'^^ a = —1/2 (i.e. 
reflected ions returning to the shock). This is the result given previously (McCIements et 
al. 1997; Laming 2001). The electromagnetic correction however has the effect of forc- 
ing the minimum wavevector to be A; = cupe/c. We require x^ > which gives k^ > 
—aujpjc?/ (a + l). As k±Ve/fte increases from the minimum wavevector decreases and 
reflected ions leaving the shock can excite the necessary waves. For k±Ve/fte ~ 1, so that 
7 ~ 0.5 and {k±/2) dl/dk± ~ -0.25 with cOpe » kv^ 

which requires ^\f,(EM)l^'^^ — 1 < a < 1/2. Figure 7 shows x^ against k±Ve/^e for « = 1/2 
and k"^ — ^^el^ (bold curve), taking iOp^. = 4 x 10^, Upi — lx 10^, and u — 560 corresponding 
to Ue — 50 and rij = 6.25 cm~^ and a magnetic held of 1.9 mG. The other curves show the 
loci of x^ for A;^ multiphed by factors of 1/2, 2, and 4 (as labeled) from this value. Apart 
from the fact that it is now ions leaving the shock that excite the waves, the value of x'^ is 
still essentially the same as that found in the cold plasma case, i.e. x'^ ~ cUpJujp^. 

Now if A; ~ ujpe/c, the approximation a; >> y/2kvi may be no longer valid in regions of 
Cas A, and we investigate the opposite limit of the plasma dispersion function, uj « y/2kvi. 
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In between these two hmits, where u; ~ y/2kvi Landau damping by thermal ions will prevent 
any wave growth. In the opposite limit (f) — > cu'^/k'^vf + .... and 

a[l + u^l/k^vl + (1 - 7) u^lJk^l] - ikj2I) dl/dk^ (1 + u^l/k^l + u^lJk^l) 



= 1- 



>e(_fc'M)/ 

(A13) 

which can be rearranged to give 



a 



(l-x^) 



,2 ,,2 



+ 1 + co^^Jk^vf + (!-/) 27 1 + co^Jk^vf + (1 - /) ^/A;^^! 

(A14) 

Again in the limit of low electron temperatures for < < 1 

k^C^-rul l-rulJQl-rulJk^vj ' ^ > 

2 

For uj^^ >> this reduces to a = k^c^^+<x>'^ ■'■ satisfies «<— l/2forA;> cUpg/c. 

In the opposite limit of high electron temperature, / — > 0, k±dl/dkj_ — > —I, and a — 
— 1/ (l + uj^/ k'^vf) giving the more usual value kmin — ^^pi/vi (see Laming 2001). However 
in this limit the lower hybrid waves do not exist, and so the appropriate wavevector for lower 
hybrid wave growth is A; = ujp^/c, arising from the electromagnetic correction. In Figure 7 
we show the generalization of this result to arbitrary electron temperature. The loci of a 
for = are plotted against k±Ve/^e, taking the following plasma parameters for Cas A, 
ujpe = 1.8 X 10"'"', ujpi = 3 X 10^, and u = 300, corresponding to Ug = 10 and rii = 1.25 cm~^, 
and a magnetic field of 1 mG. The bold curve shows the value of a for k = Wpe/c, while the 
curves show a{k±Ve/^e) for values of k successively factors of 2 greater. The dotted hue 
shows a — —1/2, remembering that necessarily a < —1/2 for ions reflected back through 
the upstream plasma with relative velocity approximately twice the shock velocity. In this 
limit wave propagating at a variety of small cosines, x to the magnetic field vector can stay 
in contact with shock reflected ions. The wave direction that dominates will be that where 
the net growth rate is maximized, as determined in section 2. 
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B. Appendix: Collision Processes 
B.l. Electron-Ion Coulomb Equilibration 

Spitzer (1978) gives the timescale for an electron distribution to relax to a Maxwellian 

3my^ {kBTef^ , . 

'^^^'^'^^ 47rV2^.eMnA ^^^^ 
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where A is the so-called plasma parameter, the ratio of largest to smallest impact parameters 
for collisions. In supernova remnants InA ~ 40. The equilibration time for ions t^g = 
teq (e, e) ^y'mi/me/Zf, and that for electron-ion equilibration is teq (e, i) — t^q (e, e) mi/me/ Z'l 
where Zi is the ion charge. Accordingly we write 

= -0.13Z=n,^ (B2) 

dt at!'^ 

which is equation (2), with AT — — Tg. We consider a fully ionized gas with Ug — Zrii 
and 

-(n.T. + neTe)=n.— + ne— = 0. (B3) 
Solving these equations yields 

dTp _ Z'^Tip Tj — XL . 
—1 = 0.13- — - (B4 
dt Z+1 ATe' 

dTj Z^Tifs Tj — Tf> . 

—1 = -0.13- — " , (B5 
dt Z + 1 ATe' 

In deriving equations (14) and (15) these expressions are averaged over the ion charge states 
in the plasma, and the expression for dT^/dt is modified by the inclusion of terms accounting 
for the change in electron density due to ionization, — (Te/rZe) {dne/dt)^^, and radiative and 
ionization losses, — (2/3neA;B) dQ/dt. Recombinations, which reduce the electron density do 
not result in an increase in the electron temperature in low density plasmas, since the energy 
of the recombined electron is radiated away, rather than being shared with the other plasma 
electrons as would be the case for three-body recombination in dense plasmas. 



B.2. Collisional Relcixation of Accelerated Electron Distribution 



The Boltzmann equation for the electron distribution function / is (see e.g. Sturrock 



1994) 



^ ^ _r A (f^ 
dt dvr \ dvr 
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(B6) 



2 dVrdVg \ dVrdVg^ 

where G = J f \v — vi \ d^vl and H = 2 f f / \v — vi \ d^vi are the Rosenbluth potentials and 
r = (47re^/mg) log (37rA) where A is the usual plasma parameter. For the collisional re- 
laxation of a suprathermal electron distribution / with a background Maxwellian fi the 
Rosenbluth potentials take the forms 
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dH 
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(B7) 
(B8) 
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where ^ {f3v) — 2/^/7r J^^ exp {—x'^) dx with = me/2kBTe. For a suprathermal electron 
distribution with >> $ (Pv) — 1. Working in one dimension then d'^G/dv'^^ — and 
dH/dv\\ — —2ne/vn. Hence 



^ = 2rnA 
dt '^dvu 




(B9) 



with / (^ii) 



imum electron ve 



- ^'Il + (^m - vf^j x^/3v^ I {vl^/2 + v^x'^/Av\;), where is the max- 
ocity in the distribution, vxi is the Alfven speed and x is the cosine of the 
angle between the magnetic field vector and k. In the cold plasma limit cos^ = ujpi/ujpe 
so vxi/x may be identified with the electron Alfven speed. The energy loss rate of the 
suprathermal electron distribution, which is equal in magnitude (but opposite in sign) to the 
heating rate of the ambient electrons, is given by 



/ 



1 2'9/ , 

-meV\\—dv\\ 

2 dt " 



Vnen'^rrie 



vl/2 + <xVK, 



2[Vm + 



log- 



Vth 



Vm + 



(BIO) 



Taking the lower limit of integration to be Vj/j = 6 x 10^ cm s~^, corresponding to an electron 
energy of 10 keV, where the ambient and accelerated electron densities are approximately 
equal, Vm — 2.5vth ~ 1.5 x 10^° cm s~^, and equation (B5) evaluates to an energy input from 
accelerated electrons to the ambient plasma of 4 x 10~^^ neu'^ ergs cm~^s~^ for vai ~ xVjn- 
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Fig. 1. — Data from the BeppoSAX MECS and PDS instruments with model bremsstrahlung 
spectra for a fully ionized oxygen plasma. The softest model spectrum is pure thermal 
bremsstrahlung. The harder spectrum has an additional 4% of electrons accelerated by 
lower hybrid waves to a maximum velocity of = 1.5 x 10^° cm s~^ (see eq. 10). The seven 
data points in bold print are taken from a rocket observation in 1968 (Gorenstein, Kellogg, 
& Gursky 1970) to illustrate the stability of at least the thermal part on the spectrum. 

Fig. 2. — Cartoon of the Gas A supernova remnant (not to scale) showing the location of 
the blast wave or forward shock, the reverse shock and the contact discontinuity between 
the reverse shocked ejecta and the forward shocked circumstellar medium. Also shown 
schematically are a few quasi-stationary fiocculi (QSFs) ahead of the forward shock and fast 
moving knots of ejecta (FMKs) ahead of the reverse shock. As the forward and reverse shocks 
encounter these density contrasts, they split into transmitted and reflected shocks. The 
reflected shocks can propagate back across the contact discontinuity. Gas A has numerous 
knots and clumps of material, making it likely that the entire shell is filled with these 
secondary shocks. 

Fig. 3. — The variation of ion temperature (dashed lines) and electron temperature (solid 
lines) for reverse shocks in Gas A in pure oxygen ejecta, against time since explosion. The 
evolution of the temperatures is plotted for ejecta encountering the reverse shock at various 
times (50, 60, 75, 100, 150, 200, and 250 years after the initial explosion), and followed to 
350 years after the explosion. The present day corresponds to i = 321 years, for an explosion 
date of 1680. The most dense ejecta is that which encounters the reverse shock earliest, and 
in this calculation undergoes radiative instability at around 120 years after explosion. 

Fig. 4. — Similar to Figure 1 but for forward shocks in Gas A at 50, 150, and 250 years 
since explosion. The elemental composition is now a mixture of N/He/H, (in mass ratios 
0.02:0.49:0.49) corresponding to that observed in the quasi-stationary fiocculi. The electron 
temperature is again the solid line, the N ion temperatures are the dashed lines, He ion 
temperatures the dotted lines, and the proton temperatures the dash-dot lines. Very little 
ion-ion or electron-ion equilibration has occurred. 

Fig. 5. — Variation of ion (dashed lines) and electron temperature (sohd lines) in pure 
oxygen ejecta reverse shocked at 50 years after the initial explosion, with electron reheating 
by lower hybrid waves commencing a further 100, 150, 200, and 250 years since reverse shock 
passage, represented by the curves labeled "b", "c", "d", and "e" respectively. Gurve "a" 
shows the initial heating and cooling following reverse shock passage, and is the same as in 
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Figure 3. The electron temperature is assumed to go no lower than 10^ K, being maintained 
at this level by residual photoionization and conduction. The electron heating is taken to 
be the same as the present day value, per unit volume. The presently observed electron 
temperature is ~ 4 x 10'' K, which is consistent with the curve with reheating commencing 
between 250-300 years since explosion, or 200-250 years since reverse shock passage. 

Fig. 6. — Data points from (Gorenstein, Kellogg, & Gursky 1970) plotted against pure 
thermal bremsstrahlung spectra corresponding to temperatures of 2 x 10^, 3 x lO'^, and 
4 X 10^ K. 

Fig. 7. — The variation of direction cosine squared, — k^k^ against k±Ve/^e for a = 1/2 
taking cOpe = 4 x 10^, Upi — 7 x 10^, and to — 560 corresponding to Ue — 50 and rii — 6.25 
cm~^ and a magnetic field of 1.9 mG. The bold curve gives the loci of x'^ for k"^ — u?^J(?. 
The other curves gives loci for equal to 1/2, 2, and 4 times this value, as labeled. For 
kxxiejVtf. < 1 direction cosines less than cupi/cupe are found, in accordance with the cold plasma 
treament. 

Fig. 8. — The variation of a with k±Ve/^e for various values of k±. The plasma parameters 
taken are cupe — 1.8 x 10^, oupi — 3x 10^, and cu — 300, corresponding to Ue — 10 and rii — 1.25 
cm~^, and a magnetic field of 1 mG. The top curve gives k±^ = uj^^jc. The successively lower 
curves are for k^ — aujpe/c where a — 1.778, 3.162, 10, and 100 respectively. 



-29- 



Table 1. Reverse Shock Parameters from Ionization Structure 



t (years) 


Q 






He (cm ^) 


Uet (cm ^s) 


approx. EM (cm ^) 


50 


0.726 


1.00e4 


1.00e4 


51.8 




5e57 


100 


0.551 


6.48e7 


1.98e7 


23.0 


2.33ell 


8e56 


150 


0.432 


5.36c8 


2.73e7 


11.8 


6.89el0 


3c56 


200 


0.348 


1.32e9 


2.37e7 


6.75 


2.20el0 


le56 


250 


0.286 


2.13e9 


1.71e7 


3.83 


6.51e9 


6c55 


300 


0.240 


2.87e9 


6.05e6 


2.00 


7.35e8 


3e55 



Table 2. Forward Shock Parameters from Ionization Structure 



t (years) 


Tn (K) 


Tne (K) 


Th{K) 




rie (cm ^) 


riet (cm ^s) 


approx. EM (cm ^) 


50 


6.30e9 


1.83e9 


4.59e8 


3.70e7 


0.96 


1.58el0 


2e55 


100 


1.02el0 


2.98e9 


7.48e8 


6.48e7 


2.73 


2.33el0 


2e56 


150 


1.13el0 


3.30e9 


8.27e8 


7.45e7 


4.58 


2.32el0 


5e56 


200 


1.05el0 


3.06e9 


7.67c8 


7.02e7 


6.20 


1.93el0 


lc57 


250 


8.99c9 


2.61e9 


6.53e8 


5.63e7 


7.44 


1.27el0 


1.5c57 


300 


7.46e9 


2.14e9 


5.37e8 


3.12e7 


8.56 


4.17e9 


2e57 



counts/cmy sec/keV 




temperature (K) 




temperature (K) 




counts/cmy sec/keV 




kiiYk^ in units (x)^l/(x)^^ 




